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ABSTRACT 

We discuss the results of N— body simulations of intermediate-mass young star clusters 
(SCs) with three different metallicities (Z = 0.01, 0.1 and 1 Z Q ), including metallicity- 
dependent stellar evolution recipes and metallicity-dependent prescriptions for stellar 
winds and remnant formation. The initial half-mass relaxation time of the simulated 
young SCs (~ 10 Myr) is comparable to the lifetime of massive stars. We show that 
mass-loss by stellar winds influences the reversal of core collapse and the expansion of 
the half-mass radius. In particular, the post-collapse re-expansion of the core is weaker 
for metal-poor SCs than for metal-rich SCs, because the former lose less mass (through 
stellar winds) than the latter. As a consequence, the half-mass radius expands faster 
in metal-poor SCs. The difference in the half-light radius between metal-poor SCs and 
metal-rich SCs is (up to a factor of two) larger than the difference in the half-mass 
radius. 

Key words: stars: binaries: general - stars: evolution - stars: mass-loss - galaxies: 
star clusters: general - methods: numerical - stars: kinematics and dynamics. 



1 INTRODUCTION 

It is well known that metallicity plays an important role 
, in the evolution of massive stars. First, it affects directly 
the luminosity and effective temperature of massive stars 
(e.g., Hurley, Pols & Tout 2000; Tumlinson & Shull 2000; 
Baraffe, Heger & Woosley 2001; Bromm, Kudritzki & Loeb 
2001). Secondly, the metallicity has a strong effect on the 
mass-loss rate by stellar winds (e.g. Kudritzki, Pauldrach & 
Puis 1987; Leitherer, Robert & Drissen 1992; Maeder 1992; 
Portinari, Chiosi & Bressan 1998; Kudritzki & Puis 2000; 
Vink, de Koter & Lamers 2001; Kudritzki 2002; Belkus, Van 
Bever & Vanbeveren 2007; Pauldrach, Vanbeveren & Hoff- 
mann 2012). This may deeply affect the evolutionary path 
in the HR diagram up to the formation of the final rem- 
nant (e.g. Heger et al. 2003a; Mapelli, Colpi & Zampieri 
2009a; Belczynski et al. 2010, hereafter B10). An interesting 
question is whether, and how much, the above effects can 
influence the overall evolution of star clusters (SCs). 

Observations indicate that there is a trend with metal- 
licity, at least for globular clusters (GCs). In fact, blue (i.e. 
generally metal-poor) GCs in the Milky Way and in some 
nearby galaxies tend to have a larger half-light radius (by 
15 — 20 per cent) than red (i.e. generally metal-rich) GCs 
(Kundu & Whitmore 1998, 2001; Kundu et al. 1999; Puzia 
et al. 1999; Larsen et al. 2001; Larsen, Forbes & Brodie 
2001; Barmby, Holland & Huchra 2002; Harris et al. 2002; 
Jordan 2004; Jordan et al. 2005, 2009; Harris 2009; Woodley 



& Gomez 2010; Strader et al. 2012). No similar studies have 
been done for young SCs (< 100 Myr) and for open clusters. 

Jordan (2004) explains the difference in the half-light 
radii of blue and red GCs as a result of mass segregation 
combined with the dependence of main sequence (MS) life- 
time on metallicity, by means of multimass isotropic Michie- 
King models. According to Jordan (2004), the observed dif- 
ference in the half-light radii does not imply a difference in 
the half-mass radii. Recent Monte Carlo (Downing 2012) and 
A-body (Sippel et al. 2012; see also Hurley et al. 2004) simu- 
lations of GCs confirm the results by Jordan (2004), finding 
that the difference in the half-light radii arises from the de- 
pendence of luminosity and stellar lifetime on metallicity. 
According to Downing (2012), there may be even a differ- 
ence in the half-mass radii, but only as a consequence of dy- 
namical interactions, mainly due to the presence of massive 
stellar black holes (BHs). Finally, Schulman, Glebbeek & 
Sills (2012) ran ./V-body simulations of young intermediate- 
mass (10 3 - 10 4 M©) SCs with different metallicities. They 
find a m 10 per cent difference also in the half-mass radius, 
between metal-poor SCs and metal-rich SCs. 

In this paper, we discuss the results of A-body simula- 
tions of young intermediate-mass SCs with different metal- 
licity and different recipes for stellar winds and remnant for- 
mation. Our aim is to investigate the core collapse and post- 
collapse evolution of young intermediate-mass SCs, to bet- 
ter understand the interplay between dynamics, metallicity- 
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dependent stellar evolution and formation of stellar rem- 
nants. 



2 THE IMPACT OF THREE-BODY 
ENCOUNTERS AND STELLAR 
EVOLUTION ON CORE COLLAPSE 

The evaporation of stars from the core of a SC removes part 
of its kinetic energy (Spitzer 1987). Since a SC has negative 
heat capacity, this leads to gravothcrmal instability and to 
the collapse of the core (Binney & Tremaine 1987). Core 
collapse in SCs is reversed mostly by three-body encounters 
(i.e. close encounters between stars and binaries). In fact, 
binaries have a energy reservoir (their internal energy, e.g. 
Binney &: Tremaine 1987), which can be exchanged with 
single stars. In particular, hard binaries (i.e. binaries with 
binding energy higher than the average kinetic energy of a 
star in the SC) tend to transfer kinetic energy to single stars 
as a consequence of three-body encounters (Heggie 1975). 
The stars that receive this kinetic energy are either ejected 
from the entire SC or remain in the periphery of the SC. 
It is worth mentioning that the more massive a binary is, 
the higher its expected encounter rate (e.g. Portegies Zwart 
2004). Thus, the most massive binaries in the SC tend to 
dominate the dynamical evolution of the system (Spitzer 
1987; Hurley 2007; Aarseth 2012; Hurley & Shara 2012). 

It has long been debated whether mass-loss by stellar 
winds and/or supernovae (SNe) is efficient in affecting core 
collapse (e.g. Angeletti & Giannone 1977, 1980; Applegate 
1986; Chernoff & Shapiro 1987; Chernoff & Weinberg 1990; 
Hurley et al. 2004; Schulman et al. 2012; Downing 2012; Sip- 
pel et al. 2012). In fact, stellar winds and SNe eject mass 
from a SC, making the central potential well shallower and 
quenching the onset of gravothermal instability. SCs with a 
broad mass-range initial mass function (IMF) undergo core 
collapse on a timescale t cc ~ 0.2 th , where th is the half- 
mass relaxation timescale (e.g. Portegies Zwart & McMil- 
lan 2002). For most GCs, th ~ 1 Gyr, whereas in young 
dense SCs th ~ 10 — 100 Myr. This means that core collapse 
in GCs is expected to occur on a timescale (much) longer 
than the lifetime of massive ( £ 20 M ) stars. Thus, the 
stages of core collapse and post-core collapse are expected 
to be barely affected by SNe and stellar winds. Instead, the 
timescale for core collapse in young dense SCs is expected 
to be of the same order of magnitude as the lifetime of mas- 
sive stars. Thus, mass-loss by stellar winds and SNe peaks 
during the epochs of core collapse and post-core collapse. 
Actually, mass-loss by stellar evolution is expected to delay 
the core collapse (quenching the gravothermal instability) 
and/or to reverse more rapidly the core collapse, depending 
on the interplay between core collapse timescale and massive 
star lifetime. 

Two further ingredients of this scenario are the depen- 
dence of mass- loss on stellar metallicity and the formation of 
stellar remnants. Stellar winds are suppressed at low metal- 
licity (e.g. Kudritzki, Pauldrach & Puis 1987; Maeder 1992; 
Kudritzki & Puis 2000; Vink, de Koter & Lamers 2001). 
Thus, metal-rich SCs are expected to lose more mass by 
stellar winds than metal-poor ones. 

Massive stars that end their life with mass higher than 
?a 40 M are expected to collapse directly into BHs, with 



no or faint SN explosion (e.g. Fryer 1999; Fryer & Kalogera 
2001). Massive metal-poor stars lose less mass by stellar 
winds, and thus are more likely to collapse directly into BHs. 
This mechanism allows to form BHs with mass higher than 
25 M (e.g. Mapelli et al. 2009a; B10). If retained inside 
the SC, these BHs become the most massive objects in the 
SC after a few tens Myr, dominating the energy budget of 
three-body encounters. 



3 METHOD 

The simulations were done using the STARLAB 1 public soft- 
ware environment (Portegies Zwart et al. 2001; see also 
Portegies Zwart & Verbunt 1996; Nelemans et al. 2001; An- 
ders et al. 2009), which allows to integrate the dynamical 
evolution of a SC, resolving binaries and three-body encoun- 
ters. In particular, we used the modified version of star- 
lab described in Mapelli et al. (2013, hereafter M13). This 
version of STARLAB includes recipes for the metallicity de- 
pendence of stellar radius, temperature and luminosity, us- 
ing the polynomial fitting formulae by Hurley, Pols & Tout 
(2000). It also includes new recipes for mass-loss by winds 
for MS stars, based on the metallicity-dependent fitting for- 
mulae given by Vink et al. (2001; see also B10). 

We added an approximate treatment for luminous blue 
variable (LBV) and for Wolf-Rayet (WR) stars. In par- 
ticular, we assume that a post-MS star becomes a LBV 
when its luminosity L and radius R satisfy the requirement 
that L/L > 6 x 10 5 and 10~ 5 (R/R e ) (L/L ) 5 > 1.0 
(Humphreys & Davidson 1994). The mass- loss rate by stel- 
lar winds for a LBV is then calculated as M = /lbv x 10~ 4 
M yr-\ where / LB v = 1.5 (B10). 

Naked helium giants coming from stars with zero age 
MS (ZAMS) mass m Z AMS > 25 M (e.g., van der Hucht 
1991 and references therein) are labelled as WR stars in the 
new version of the code and undergo a mass-loss rate by 
stellar winds defined by M = 1O" 13 (L/L ) 1 ' 5 (Z/Zq) 13 M 
yr _1 , where f3 = 0.86. This formula was first used by B10, 
and is a combination of the Hamann & Koesterke (1998) 
wind rate estimate (taking into account WR wind clump- 
ing) and Vink & de Koter (2005) wind Z-dependence for 
WR stars. Stellar winds in asymptotic giant branch (AGB) 
stars are modelled as in the standard version of STARLAB 
(Portegies Zwart & Verbunt 1996), i.e. the code does not 
include any recipes for metallicity-dependent stellar winds 
in AGB stars. 

The formation of stellar remnants is implemented as de- 
scribed in M13. In particular, BH masses for various metal- 
licities follow the distribution described in fig. 1 of M13 (see 
also Fryer et al. 1999; Fryer & Kalogera 2001; B10; Fryer et 
al. 2012). If the final mass 2 of the progenitor star is > 40 
M , we assume that the SN fails and that the star collapses 
quietly to a BH. The requirement that man > 40 M im- 
plies that only stars with ZAMS mass > 80 and > 100 M , 
can undergo a failed SN at Z = 0.01 and 0.1 Z , respec- 
tively. If mfln > 40 M , the mass of the BH is derived as 
msH = mco + Zcoii (mne + mn), where mco is the final mass 

1 http : //www. sns . ias . edu/~ starlab/ 

2 We call 'final mass', m.fl n , of a star the mass bound to the star 
immediately before the collapse. 
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Table 1. SC properties in initial conditions. 



Parameter 


Values 


Wo 


5 


TV* 


5500 


r c (pc) 


0.4 


c 


1.03 


IMF 


Kroupa (2001) 


m-min (Mq) 


0.1 


"imax (Mq) 


150 


fpB 


0.0, 0.1 


Z(Z ) 


0.01, 0.1, 1.0 



Wo- central dimcnsionlcss potential in the King (1966) model; 
TV*: number of stars per SC; r c : initial core radius; 
c = log io (ft /V c ): concentration (rt is the initial tidal radius); 
fn m i n and m max : minimum and maximum simulated stellar 
mass, respectively; /pg: fraction of primordial binaries, defined 
as the number of primordial binaries in each SC divided by the 
number of 'centres of mass' (CMs) in the SC. In each simulated 
SC, there are initially 5000 CMs, among which 500 arc 
designated as 'binaries' and 4500 are 'single stars' (see Downing 
ct al. 2010 for a description of this formalism). Thus, 1000 stars 
per SC are initially in binaries. 

of the Carbon Oxygen (CO) content of the progenitor, while 
mHc and mH are the residual mass of Helium (He) and of 
Hydrogen (H), respectively. / co ii is the fraction of He and H 
mass that collapses to the BH in the failed SN scenario. The 
value of /con is uncertain. We assume / co ii = 2/3 to match 
the maximum values of mrm at low Z derived by B10. In this 
scenario, BHs with mass up to ~ 80 Mq (~ 40 Mq) can form 
if the metallicity of the progenitor is Z ~ 0.01 Z (Z ~ 0.1 
Z Q ). BHs that form from quiet collapse are assumed to re- 
ceive no natal kick (see Fryer et al. 2012). For BHs that 
form from a SN explosion, the natal kicks were drawn from 
the same distribution as neutron stars but scaled with the 
square root of the mass (see M13 for details). 

We assume that the mass lost by stellar winds and SNe 
is immediately removed from the simulation. This assump- 
tion is correct for SN ejecta and also for the winds of massive 
stars, which are expected to move fast ( £ 2000 km s _1 for 
the O stars, e.g. Muijres et al. 2012; > 1000 km s _1 for 
the WR stars, e.g. Vink & De Koter 2005; Martins et al. 
2008) with respect to the escape velocity of the simulated 
SCs ( £ 10 km s -1 ). Stellar winds by AGB stars have much 
smaller velocities (w 10 — 20 km s~\ Loup et al. 1993), but 
still sufficiently high to escape from our simulated SCs. Fur- 
thermore, we show in Section 4 that AGB stars do not play 
an important role for the results presented in this paper. 

3.1 Initial conditions and simulation grid 

The main properties of the simulated SCs are the same as de- 
scribed in M13. In particular, we focus on intermediate-mass 
(Mtot = 3000-4000 Mq) young (< 100 Myr) SCs. We as- 
sume a spherical King profile with central dimensionless po- 
tential Wo — 5 (King 1966), initial core radius r c = 0.4 pc, 
concentration c = log 10 (r t /r c ) = 1.03 (where r t is the tidal 
radius). The resulting half- mass radius is rh m ~ 0.8 — 0.9 
pc. The basic SC properties are listed in Table 1. 



The initial centres of mass (CMs) of the particles in 
each simulation arc 5000. Each CM corresponds either to 
a single star or to the CM of a binary. This formalism is 
commonly used in simulations of SCs including primordial 
binaries (e.g., Downing et al. 2010). Note that the fraction 
of primordial binaries /pb is defined as the number of bina- 
ries, divided by the total number of CMs. Thus, /pb =0.1 
means that there are 500 binaries over 5000 CMs (i.e. 5500 
stars). The single stars and the primary members of a binary 
follow a Kroupa initial mass function (IMF, Kroupa 2001), 
with minimum and maximum mass equal to 0.1 and 150 
Mq, respectively. The masses of the secondary stars (7712) 
are generated according to a uniform distribution between 
0.1 mi and mi (where mi is the mass of the primary). The 
initial semi-major axis a of a binary is chosen from a distri- 
bution /(a) tx \ja (Sigurdsson & Phinney 1993; Portegies 
Zwart &; Verbunt 1996), consistent with the observations of 
binary stars in the Solar neighbourhood (e.g. Kraicheva et 
al. 1978; Duquennoy & Mayor 1991). We generate a between 
R and 10 5 Rq , but discarding systems where the distance 
between the two stars at the pericentre is smaller than the 
sum of their radii (Portegies Zwart, McMillan & Makino 
2007). The initial eccentricity e of a binary is chosen from a 
thermal distribution /(e) = 2e, in the — 1 range (Heggie 
1975). 

The half-mass relaxation time for the simulated SCs is 
t h ~ 10Myr(r h /0.8pc) 3/2 (Mtot/3500M q ) 1/2 . Thus, the 
core collapse time (Portegies Zwart & McMillan 2002) is 
t cc « 2 — 3 Myr (th/10 Myr). We integrate the evolution of 
the SCs for the first 100 Myr, i.e. the epoch when the in- 
terplay between strong dynamical interactions and massive 
stellar evolution is more important. The properties of the 
grid of simulations are summarized in Table 2. For each SC 
model, we run a number TV rc of single realizations (chang- 
ing only the random seeds), to filter out the fluctuations 
associated with each single realization and to get an 'aver- 
age' model. Runs Al, B and C are our fiducial runs (i.e. 
the models described in M13) and differ only for the metal- 
licity (Z — 0.01, 0.1 and 1 Z , respectively 3 ). The runs 
labelled with Ai (where i = 2,3, 4, 5) have the same metal- 
licity (Z — 0.01 Zq) as our fiducial runs Al, and differ from 
Al for other parameters. In particular, (i) in A2 /pb = 0, 
(ii) in A3 the maximum allowed BH mass is m.BH, max = 25 
Mq, (iii) in A4 all BHs receive a natal kick «kick = 10 3 km 
s _1 (set to an unrealistically high value to eject all BHs from 
the SC, for comparison with Downing 2012), and (iv) in A5 
no mass-loss by stellar winds was implemented. 

The properties of the simulated SCs (total mass, num- 
ber of stars, core density, core and half-mass radius) are con- 
sistent with the properties of observed young intermediate- 
mass SCs (see e.g. the recent review by Portegies, McMillan 
& Gieles 2010; see also Hillenbrand & Hartmann 1998; Dias 
et al. 2002; Portegies Zwart 2004; Pfalzner 2009; Kuhn et al. 
2012). Finally, our simulations do not include recipes for the 
tidal field of the host galaxy. Accounting for the tidal field 
may increase the fraction of mass lost and even transform 
the SCs into unbound associations (e.g. Gieles & Portegies 
Zwart 2011). The effect of tidal fields will be added and 
discussed in forthcoming papers. In this paper, we decided 



In our simulations, we assume Z = 0.019. 
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Table 2. Simulation grid. 



Name 




Z 


Stellar winds 


BH kick 


/PB 








(Zq) 








(M Q ) 


Al 


100 


0.01 


YES 


LOW 


0.1 


80 


B 


100 


0.1 


YES 


LOW 


0.1 


40 


C 


100 


1.0 


YES 


LOW 


0.1 


23 


A2 


100 


0.01 


YES 


LOW 


0.0 


80 


A3 


50 


0.01 


YES 


LOW 


0.1 


25 


A4 


22 


0.01 


YES 


HIGH 


0.1 


80 


A5 


22 


0.01 


NO 


LOW 


0.1 


110 



Column 1: Name of the set of runs; Column 2, N [c : number of 
random realizations per each SC model; Column 3: metallicity; 
Column 4: YES/NO distinguishes between models in which 
stellar winds are/are not included; Column 5: LOW/HIGH 
distinguishes between models where BH natal kicks are 
described as in M13 (i.e., no kick is assigned to BHs with mass 
> 25 Mq; whereas the natal kicks were drawn from the same 
distribution as neutron stars but scaled with the square root of 
the mass for BH masses < 25 Mq), and models where a natal 
kick velocity = 10 3 km s -1 was assigned to all BHs (this 

velocity was purposely set to an unrealistically high value, to 
eject all BHs from the SC). Column 6, fpB- fraction of 
primordial binaries; Column 7: »tibh, max is the maximum 
possible mass for BHs. Runs Al, B, C, A2 and A3 were already 
presented in M13. 
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Figure 1. Number of core-collapse SNe in our simulations as a 
function of time for Z = 0.01 Zq (cross-hatched red histogram, 
runs Al), Z = 0.1 Zq (black empty histogram, runs B) and Z = 1 
Zq (hatched blue histogram, runs C). Vertical dashed line: time 
below which SNe are failed for Z = 0.01 Z and Z = 0.1 Zq, in 
our simulations. In this and in all the figures of this paper, colours 
are available in the online version. 



not to include tidal fields because we want to look at the 
intrinsic properties of the simulated SCs. 



4 RESULTS 

The simulated SCs have initial half-mass relaxation time 
th ~ 10 Myr and core collapse time t cc « 3 Myr (see the 
previous section). These timescales are of the same order of 
magnitude as the lifetime of the most massive stars. The 
lifetime of a 30 M star is ~ 6 Myr, and stellar winds are 
relatively inefficient for smaller stellar masses (excluding the 
AGB phase). Thus, the most intense phase of mass-loss by 
stellar winds coincides with the collapse and re-expansion of 
the SC core: thus, the effects of stellar winds on the struc- 
tural properties of the SC should be maximal in our simu- 
lations. 

Core-collapse SNe span a longer time-range with respect 
to stellar winds. Fig. 1 shows that the rate of SNe is max- 
imum at t ~ 3 — 10 Myr, remains quite high up to w 50 
Myr, and drops at t > 50 Myr. The few SNe at t > 50 Myr 
involve blue straggler stars (BSSs), which are the result of 
either mass transfer or a merger between two MS stars (e.g. 
Mapelli et al. 2004, 2006, 2009b). For 50 < f/Myr < 100, 
the number of SN explosions involving BSSs is ~ 0.35 per 
SC (regardless of the metallicity), corresponding to a mass- 
loss of ~ 2.5 Mq per SC. Since the total mass-loss in the 
time interval between 50 and 100 Myr is « 3.5 per cent of 
the initial total mass of the SC (Mtot), i.e. ~ 100 - 140 
M per SC, the mass lost through SN explosions of BSSs is 
negligible. Thus, after t ^ 50 Myr, the evolution of the core 
is dominated by three-body interactions and stellar winds 
of AGB stars. 




t (Myr) 



Figure 2. Cumulative mass-loss by stellar winds and SNe 
(Mlost) normalized to the initial total mass of the SC (Mtot) 
as a function of time for the first 6.5 Myr. In the insert: 
Mlost/Mtot is shown for 100 Myr. Solid red line: Al (Z = 
0.01 Zq); dashed black line: B (Z = 0.1 Zq); dotted blue line: C 
(Z = IZq). Each line in this Figure is the median value of 100 
simulated SCs. 
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Figure 3. Core radius (r c ) as a function of time for the three 
considered metallicities: Z = 0.01 Zq (solid red line, Al), Z = 
0.1 Zq (dashed black line, B) and Z = 1 Zq (dotted blue line, C). 
Each line in this Figure is the median value of 100 simulated SCs. 
In the insert: zoom of the first 7.5 Myr. 

For comparison, most GCs have half-mass relaxation 
timescales and core collapse timescales that are a factor of 
£ 10 longer (e.g. Portegies Zwart 2004; Portegies Zwart, 
McMillan & Gieles 2010), indicating that the peak of stellar 
wind activity and that of SN explosions ended well before 
the beginning of the core instability phase. 

4.1 Metallicity dependence 

Fig. 2 shows the cumulative mass lost by the SC (because 
of both stellar winds and SNe), as a function of time, for 
runs Al, B and C. The effect of different metallicities is 
maximum between 1 and 4 Myr, when the mass lost at Z — 1 
Zq is up to 100 and 10 times higher than the mass lost at 
Z — 0.01 Zq and Z = 0.1 Zq, respectively. The mass lost in 
these early phases of the SC life is a relatively small fraction 
(< 0.1) of the total SC mass, but is large when compared to 
the core mass (which is ^0.1 Mtot). As the most massive 
stars already sank to the core through dynamical friction 
by the time of core collapse, the effect of early mass-loss is 
particularly strong in the SC core. For example, at t = 3.3 
Myr the mass-loss is ~ 30 per cent of the core mass at Z — 1 
Zq, and only ~ 1 per cent at Z — 0.01 Zq. 

Fig. 3 shows the behaviour of the core radius r c as a 
function of time for runs Al, B and C (i.e. for the three con- 
sidered metallicities). The collapse is so fast (t cc ~ 3 Myr) 
that it occurs almost at the same time for all the consid- 
ered metallicities. The effect of metallicity in Fig. 3 appears 
immediately after the collapse, during the first phase of re- 
expansion: the core radius in metal-rich SCs expands more 
than in metal-poor SCs, as mass-loss by stellar winds in 
metal-rich SCs removes more matter from the core poten- 
tial well. At t w 6 Myr, the mass-loss by stellar winds is 
nearly over, but the differences among core radii at different 
metallicity remain almost constant up to t ~ 30 Myr. 



Figure 4. Half-mass radius (rhm) as a function of time for the 
three considered metallicities: Z = 0.01 Zq (solid red line, Al), 
Z = 0.1 Zq (dashed black line, B) and Z = 1 Zq (dotted blue 
line, C). Each of the aforementioned three lines in this Figure is 
the median value of 100 simulated SCs. The dot-dashed magenta 
line shows the analytical prediction. In the insert: zoom of the 
first 7.5 Myr. 

The contribution of SN explosions lasts for a longer time 
(t ~ 50 Myr, Fig. 1). In our models, the mass- loss by SN 
explosions does not depend on metallicity. The only impor- 
tant effect of metallicity on SNe is that failed SNe can take 
place only at low metallicity ( £ 0.1 Zq). All the failed SNe 
occur at t < 3.5 Myr (dashed line in Fig. 1), as only the 
most massive stars (see Section 3) can undergo a failed SN. 
Thus, the occurrence of failed SNe in low-metallicity SCs 
enhances the difference between metal-rich and metal-poor 
SCs, in the early phase of core-collapse reversal. 

At times t > 50 Myr, even the mass-loss by SNe is 
over: the only process that may affect significantly the later 
evolution of the core radius is represented by three-body 
encounters. The differences between metal-rich and metal- 
poor SCs tend to be quenched at late times. 

The half-mass radius rhm (Fig- 4) remains almost con- 
stant during core collapse (as expected, e.g. Elson, Hut & 
Inagaki 1987), while it starts increasing after the reverse of 
the core collapse. The behaviour of rhm {rhm oc t 2//3 ) is in 
agreement with simple analytical predictions (e.g. Elson et 
al. 1987). The post-core collapse value of rhm for metal-poor 
SCs is systematically larger than that for metal-rich SCs, in 
agreement with Schulman et al. (2012). The reason is that 
the reversal of core-collapse is slower for metal-poor SCs. 
This implies that metal-poor SCs maintain a higher core 
density in the late core-collapse phase and in the early core- 
collapse reversal (see Fig. 5 for the evolution of core stel- 
lar density 4 ). Since the rate of three-body encounters scales 
approximately with the stellar mass density (e.g. Sigurds- 
son & Phinney 1993), metal-poor SCs have a higher rate 

4 The core mass (number) density of stars was approximated as 
the total mass (number) of stars in the core divided by rj? . 
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Figure 5. Top panel: core mass density of stars (p) as a function 
of time for the three considered metallicities: Z = 0.01 Zq (solid 
red line, Al), Z = 0.1 Z (dashed black line, B) and Z = 1 Z Q 
(dotted blue line, C). Each line in this Figure is the median value 
of 100 simulated SCs and is normalized to poi i- e - the core mass 
density of stars in the case of Z = 0.01 Zq. Bottom panel: core 
number density of stars (n) as a function of time for the three 
considered metallicities: Z = 0.01 Zq (solid red line, Al), Z = 
0.1 Zq (dashed black line, B) and Z = 1 Z (dotted blue line, C). 
Each line in this Figure is the median value of 100 simulated SCs 
and is normalized to no, i.e. the core number density of stars in 
the case of Z = 0.01 Zq. 



Figure 6. Top panel: half-light radius (r^i) as a function of time 
for the three considered metallicities: Z = 0.01 Zq (solid red line, 
Al), Z = 0.1 Zq (dashed black line, B) and Z = 1 Zq (dotted 
blue line, C). Each line in this Figure is the median value of 100 
simulated SCs, and is normalized to the median half-light radius 
of SCs with Z = 0.01 Zq (^mo)- Bottom panel: half-mass radius 
(?"lim) & s a function of time for the three considered metallicities: 
Z = 0.01 Zq (solid red line, A), Z = 0.1 Zq (dashed black line, 
B) and Z = 1 Zq (dotted blue line, C). Each line in this Figure is 
the median value of 100 simulated SCs, and is normalized to the 
median half-mass radius of SCs with Z = 0.01 Zq (r^mo)- 



of three-body encounters than metal-rich SCs. Three-body 
encounters pump kinetic energy in the SC halo (mainly in 
the form of stars ejected in the outskirts of the SC), and are 
responsible for the expansion of rhm (e.g., Elson et al. 1987). 

The difference at t £ 40 Myr between the averaged rh m 
of Z = 0.01 Zq SCs and that of Z = 1 Z© SCs is ~ 10 per 
cent, surprisingly similar to the difference observed between 
red and blue GCs. On the other hand, we stress that the evo- 
lution of our simulated young SCs is very different from the 
evolution of GCs. Fig. 6 shows that the difference among 
the half-light radii at different metallicity is larger (w 20 per 
cent, although with larger fluctuations) than the difference 
among the half-mass radii. This results from the combina- 
tion between mass segregation, and the metallicity depen- 
dence of the adopted stellar luminosity function. In our sim- 
ulated SCs, mass segregation is very efficient, as shown by 
the fact that the core collapse occurs on a timescale shorter 
than the half-mass relaxation timescale (i.e., it is driven by 
dynamical friction, e.g. Portegies Zwart 2004). Thus, the re- 
gion inside the core radius is dominated by massive stars 
and remnants, while most low-mass stars are in the outer 
regions. 

According to the fitting formulae by Hurley et al. 
(2000), a solar metallicity MS star with mass ^$ 15 Mq is 
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Figure 7. ZAMS luminosity (-Lzams) versus ZAMS mass for 
Z = 0.01 Zq (solid red line), Z = 0.1 Z (dashed black line) and 
Z = 1 Z (dotted blue line). LzAMS i s normalized to Lzams(-^ = 
0.01 Z Q ), i.e. to the ZAMS luminosity for Z = 0.01 Z Q . 
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Figure 8. Core radius (r c ) as a function of time for Z = 0.01 Zq. 
Solid red line: Al, fiducial model with /pb = 0.1 (the same as in 
Fig. 3); long dash-dotted green line: A2, the same as the fiducial 
model, but with /pb = 0.0; short dash-dotted magenta line: A3, 
m BH max = 25 Mq; dashed black line: A4, BHs ejected by natal 
kick; dotted blue line: A5, no mass-loss by stellar winds. Each line 
in this Figure is the median value of a number of simulated SCs 
(see Table 2). In the insert: zoom of the first 7.5 Myr. 

fainter than a sub-solar metallicity MS star with the same 
mass. This is apparent from Fig. 7, where we compare the 
ZAMS luminosity (Lzams) for stars with different metal- 
licity. A similar difference persists during the entire MS. 
Because of this difference in the luminosity function, and 
because our simulated SCs are mass-segregated, the light 
distribution tends to be more concentrated in metal-rich SCs 
than in metal-poor SCs. 

4.2 Other effects related to stellar evolution 

In this section, we estimate the impact on core collapse and 
post-core collapse phases of other effects connected with stel- 
lar evolution, including stellar remnants. To maximize the 
contribution of BHs to the reverse of core collapse, we con- 
sider SCs with Z — 0.01 Zq, where the BH mass can be as 
high as ~ 80 M . 

The comparison between runs Al (with /pb = 0.1) and 
A2 (with /pb = 0) shows that the existence of primordial 
binaries has almost no effect on both the core radius (Fig. 8) 
and the half-mass radius (Fig. 9). In fact, the very hard 
binaries needed to reverse the core collapse form by 3-body 
capture during the collapse phase even in the /pb = runs. 

The difference between runs Al and A3 is the maximum 
mass of BHs: 80 and 25 Mq for runs Al and A3, respectively. 
This difference was achieved by ejecting impulsively a large 
fraction of the stellar mass at the time of BH formation in 
runs A3. Thus, the core radius in runs A3 increases much 
faster than that in runs Al because of this sudden mass-loss 
(Fig 8). On the other hand, the half-mass radius (Fig 9) in- 
creases initially faster in runs A3 than in runs Al, but then 
slows down, because of the lower core density and because 
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Figure 9. Half-mass radius (rjun) as a function of time for Z = 
0.01 Zq. Solid red line: Al (the same as Fig. 4); long dash-dotted 
green line: A2; short dash-dotted magenta line: A3; dashed black 
line: A4; dotted blue line: A5. Each line in this Figure is the 
median value of a number of simulated SCs (see Table 2). In the 
insert: zoom of the first 7.5 Myr. 

of the absence of massive BHs. In fact, massive BH binaries 
have a larger cross-section for three-body encounters, have 
(on average) a larger reservoir of internal energy and, thus, 
are more efficient in ejecting stars after three-body encoun- 
ters (e.g. Downing 2012). 

A similar but stronger trend can be observed by com- 
paring runs Al and A4. In runs A4 all BHs are removed from 
the simulation at birth, by assigning to them a natal kick of 
10 3 km s" 1 . Immediately after core collapse, the core radius 
(Fig 8) and the half-mass radius (Fig 9) in runs A4 become 
significantly larger than in runs Al. In fact, the mass of all 
the BHs is completely lost from the potential well of the SC. 
On the other hand, the removal of all the BHs implies that 
less energy can be exchanged via three-body encounters. For 
this reason, the increase of the half-mass radius slows down 
at t £ 40 Myr. 

Finally, runs A5 represent the most extreme case, where 
mass-loss by stellar winds is switched off and the maximum 
BH mass is 110 Mq (the only difference between the BH 
mass and the ZAMS mass of the progenitor star comes from 
the recipes for the direct collapse of the star into a BH, see 
Section 3). If mass- loss is switched off, the core collapse is 
much more dramatic: it lasts for 20 Myr before that three- 
body encounters can reverse it (Fig 8). The half-mass radius 
(Fig 9) expands dramatically after core collapse, because of 
the amount of energy pumped into the halo by three-body 
encounters. 



5 CONCLUSIONS 

In this paper, we discussed the effects of metallicity on 
the core collapse and post-core collapse phase in young 
intermediate-mass SCs. This was done by means of iV-body 
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simulations including recipes for metallicity-dependent stel- 
lar evolution, stellar winds and formation of stellar rem- 
nants. 

In the simulated SCs, core collapse is almost coincident 
with the peak of mass-loss by stellar winds. Metal-rich SCs 
lose more mass by stellar winds than metal-poor SCs: the 
reversal of core collapse is faster and stronger in the former 
with respect to the latter. On the other hand, the difference 
in the core radius among metal-rich and metal-poor SCs 
decreases as soon as mass losses by massive stars are over. 
In the later stages of SC life, the core evolution is ruled 
mainly by three-body encounters. 

Since the reversal of core collapse is slower in metal- 
poor SCs, the half-mass radius expands more in metal-poor 
SCs than in metal-rich SCs. The maximum difference among 
the half-mass radii of Z — 0.01 and Z = 1 Zq SCs is ~ 
10 per cent. When considering the half-light radii rather 
than the half-mass radii, the difference at late times is larger 
(« 20 per cent). This is a consequence of the metallicity 
dependence of the adopted stellar luminosity function. 

We also checked the effects of other aspects of stellar 
evolution. When stellar winds are completely suppressed, 
the core collapse phase lasts much longer and the half-mass 
radius in the post-core collapse phase can be even 20 per cent 
larger than in the fiducial model for the same metallicity. 
If all the BHs are ejected by natal kick, the core initially 
expands faster because of the impulsive mass-loss, but all 
the massive binaries are lost from the SC, quenching the 
effects of three-body encounters. 

In the last few years, young SCs were at the centre 
of important observational campaigns (e.g. Bica et al. 2003; 
Mercer et al. 2005; Borissova et al. 2011; Bianchi et al. 2012; 
Richards et al. 2012; see Portegies Zwart et al. 2010 for a 
recent review), and are one of the main targets of the on- 
going Gaia ESO survey (Gilmore et al. 2012). Thus, com- 
paring the predictions of our simulations with the current 
and forthcoming data about metallicity and half-light radii 
of young SCs will give new insights on the formation and 
dynamical evolution of SCs. 
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